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We investigate the dynamic behavior of a Bose-condensed gas of alkah atoms interact- 
ing with repulsive forces and confined in a magnetic trap at zero temperature. Using the 
Thomas-Fermi approximation, we rewrite the Gross-Pitaevskii equation in the form of the 
hydrodynamic equations of superfluids. We present solutions describing large amplitude os- 
cillations of the atomic cloud as well as the expansion of the gas after switching off the trap. 
We compare our theoretical predictions with the recent experimental data obtained at Jila 
and MIT. 
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The dynamic behavior of Bose-condensed gases of al- 
kali atoms in magnetic traps has been the object of recent 
experiments at Jila and MIT [|j|J|l. Absorption 

images of the atomic cloud provide quantitative informa- 
tion on the dynamics of the expansion after switching 
off the trap as well as accurate data for the frequencies 
of the collective excitations. The experimental results 
reveal that the role played by the interatomic forces is 
important in these systems and can not be ignored for a 
quantitative understanding of the data. This opens new 
challenging tasks for theoretical investigation. 

The natural theory to investigate the dynamic behav- 
ior of a nonuniform Bose condensate, at T = 0, is given 
by the time-dependent Gross-Pitaevskii (GP) equation 
for the condensate wavefunction \E'(r,t): 
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The coupling constant g is proportional to the s-wave 
scattering length a through g = inh^a/m,. In the fol- 
lowing we will discuss the case of repulsive interactions 
(a > 0). The anisotropic trap is represented by the 
confining potential Vext, which is chosen in the form 



where = x, y, z. So far 



the experimental traps have cylindrical symmetry and 
hence are characterized by the radial frequency ll!±_ = 
LOx = and the asymmetry parameter A = lOzj^L- The 
ground state configuration [p|,p|JTc|] as well as the proper- 
ties of small oscillations near equilibrium have 
been the object of systematic investigation starting from 
Eq. (Q). A few calculations in the nonlinear regime have 
been also carried out ||l|,|l|,|l|,0,|8| . 

In the present paper we discuss several features of the 
nonlinear behavior of the system, by solving Eq. (^) in 
the large N limit, where N = J'dr|5'(r, <)p is the num- 



ber of atoms. In this limit, it is possible to derive almost 
analytic results of the GP equation, thereby simplifying 
the numerical analysis and allowing for a systematic in- 
vestigation of important phenomena. These include, the 
explicit time evolution of the condensate (shape of pro- 
files, aspect ratio etc.) during the expansion and the de- 
pendence of the collective frequencies on the amplitude 
of the oscillation. 

The effective strength of the interatomic forces in the 
GP equation is fixed by the adimensional parameter 
Na/a^i, where a^i = \/'hl {mujQi) is the harmonic os- 
cillator length. When this parameter is much larger than 
1, the repulsion makes the system much wider than the 
noninteracting configuration, yielding a rather smooth 
profile. In such conditions, the equilibrium results from 
a balance between the external potential and the repul- 
sive interaction, the kinetic energy playing a minor role. 
For large values of Na/aQi, one can then neglect the ki- 
netic energy term in (|l|). This yields the Thomas-Fermi 
(TF) approximation for the ground state: 
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when fj, > yea;t(r) and poi'r) = elsewhere. The chemical 
potential /x is fixed by the normalization of the density 
to the number of particle N: 
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The Thomas-Fermi approximation (|^) works in an excel- 
lent way for the configurations realized at MIT, where N 
is of the order of 1 million atoms and more. Conversely, in 
the Jila experiments of Refs. the number of atoms 

is smaller (10^-10^) and the TF approximation provides 
only a semi-quantitative description. 
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Neglecting the kinetic energy pressure term in the time 
dependent GP equation corresponds to taking into ac- 
count the effect of the kinetic energy operator in (|l]) only 
on the phase of the order parameter ^E*. This permits to 
rewrite (uh in the useful hydrodynamic form p3[ : 



^p + V- (vp) = 

TO^V + V I Vext +9P+ — 



(4) 
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where density and velocity are defined by p = and 
V = (\['*V\E'-*V^'*)?i/(2mip). Equation is the usual 
equation of continuity, while (|^) establishes the irrota- 
tional nature of the superfluid velocity. It is immedi- 
ate to verify that the equilibrium configuration given by 
Eqs. (Hjq) (v = and dp/dt = 0) coincides with the TF 
result (^. 

The hydrodynamic equations have been recently 
shown to provide the correct frequencies of the normal 
modes of the condensate in the large TV limit . With 
respect to the full solution of the GP equation, which 
includes the effect of the kinetic energy pressure, the ap- 
proach based on (Q-H) has the major advantage of pro- 
viding an algebraic expression for the dispersion relation 
of the elementary excitations. The resulting predictions 
quite well compare with both the Jila and MIT ||] 
experiments. 

For nonlinear time dependent motions, which are the 
object of the present work, almost analytic solutions can 
be found starting from equations (|^^. In fact they ad- 
mit exact solutions of the form 

p(r, t) = a^{t)x'^ + ay{t)y^ + a,{t)z^ + ao{t) (6) 
V = ^W[a4t)x^ + ay{t)y^ + a^{t)z^] . (7) 

Equation is restricted to the region where p> and 
the coefficient uq is fixed by the normalization of the den- 
sity: ao = -{15N/8TT)'^/^{axayaz)'^^^. The time depen- 
dent coefficients and ai obey the following coupled 
differential equations: 



di + 2aiai + Ui aj = 
j 

di + af + + {2g/m)ai — , 
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with i,j — x,y, z. One can use the three equations (||) to 
express Ui in terms of hi/ai] the solution is greatly sim- 
plified by introducing the new adimensional variables bi 
defined by = —mujQ^{2gbxbybzb^)~^. With this choice, 
Eqs. (§) reduce to ai = bi/bi and Eqs. become 



bi + Lolibi - LolJ (bib^byb^) = 



(10) 



The second and third terms of (jl^) give the effect of the 
external trap and of the interatomic forces, respectively. 



It is worth noticing that, using the new variables bi, the 
equations of motion do not depend on the value of the 
coupling constant g. This is a typical feature characteriz- 
ing the large N behavior of the GP equation. The mean 
square radii and velocities of the atomic cloud can be 
easily expressed in terms of bi: 
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where /i is given by (^. 

The solutions (^[-Q) are well suited to describe both the 
oscillations around the ground state and the problem 
of the expansion of the gas after switching off the confin- 
ing potential. Equations ( [To| ) have been already derived 
by other authors |l^Jl^, using the formalism of scaling 
transformations, and applied to study several dynamic 
phenomena. In the present approach the same equations 
emerge as an exact solution of the hydrodynamic equa- 
tions of supcrfluids (|^||) and will be used to investigate 
nonlinear oscillations and the expansion in both the Jila 
and MIT traps. 

In order to apply the above formalism to study the ex- 
pansion of the cloud, let us suppose that at t = the 
system is in its equilibrium configuration. Comparing 
Eq. (||) with the TF ground state density (||) and sub- 
stituting ai in terms of bi, one finds bi = 1. One has 
also bi = at equilibrium. The external potential is then 
suddenly switched off and the system starts expanding. 
This corresponds to solving the equations with the 
second term set equal to zero: 



bi - ioli/ihb^byb^) = 



(13) 



For an initially spherical configuration the expansion 
will proceed isotropically. In the presence of anisotropy, 
the expansion (and consequently the asymptotic veloc- 
ities) will be instead faster in the direction where the 
repulsive forces (proportional to the gradient of the den- 
sity) are stronger. For an axially deformed trap with 
A — ujz/i^i_ > 1 this will occur in the axial direction, 
while if A < 1 it will occur in the radial direction. The 
ratio 



Rr{t) = ^{z^)/{x^) = X~'bAt)/b,{t) (14) 

of the radii in the two different directions is called the 
aspect ratio in co-ordinate space. One can also define 
the aspect ratio of velocities 



R,{t)^^{v^_)/{vl)=X-%{t)/b,it), (15) 

whose deviations from unity refiect the anysotropy of 
the velocity distribution. This anysotropy represents 
a crucial feature of Bose condensates. Asymptotically 
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{t — > oo) the two aspect ratios Rr and Ry converge to 
the same value R. However for finite values of the expan- 
sion time they can behave quite differently. In fact, the 
velocities in both radial and axial directions reach their 
asymptotic values very rapidly, since the time scale for 
acceleration is very short; vice versa the radii approach 
their asymptotic behavior, oc Vit, much more slowly. This 
is an important feature to take into account in the anal- 
ysis of experimental data. 

In Figs. |l| and |^ we show the results of the two as- 
pect ratios Rr and R^ obtained by solving numerically 
Eqs. ( JTgI ) for two different sets of frequencies woi, cor- 
responding to the Jila [^,|| and MIT § traps, respec- 
tively. The different behavior of the two aspect ratios 
Rr and R^ is evident in both cases. In the same figures 
we have also shown the predictions of the noninteracting 
harmonic oscillator model, which gives (vf) — hujQi/{2m) 
and (r^) = (l/2)aoj + {vf)t^. The figures point out 
very clearly the role of two-body interactions which mod- 
ify both the timescale of the expansion process and the 
asymptotic value of the aspect ratio. The comparison 
with the experimental data for Rr should be however 
taken with care. In the case of the Jila experiments, N is 
of the order of 10'^ ~ 10*, so that the TF approximation 
is expected to be rather crude; moreover, the points in 
Fig. are taken from a gaussian fit to the spatial distribu- 
tion of the atoms and they can be significantly corrected 
by using different fitting functions | p5| , p0| . In the case 
of the MIT experiments, the points in Fig. |^ correspond 
to our estimate of the aspect ratio, extracted from the 
time-of- flight images in Fig. 1 of Ref . Q . 

In Fig. ^ we report the asymptotic aspect ratio R as 
a function of A. This curve has been recently calculated 
also in Ref. ||l^. In the limit A ^ the aspect ratio 
approaches the value (7r/2)A [|l^. For comparison we 
also show the predictions of the noninteracting harmonic 
oscillator model, given by Rho = a/A- 

Another important quantity to discuss is the release 
energy Erei, defined as the energy per particle of the 
system after the switch-off of the trap. This energy is 
given by the sum of the kinetic and interaction energy 
of the atoms and it is conserved during the expansion, 
being finally converted entirely into the kinetic energy 
of the expanding cloud. In the present formalism the 
release energy is a first integral of equation ([l3| ) and can 
be written as 
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where fj, is given by (^. At the beginning of the expan- 
sion Erei coincides with the 2-body interaction energy of 
the TF ground state. The comparison between Eq. ( p^ ) 
and the numerical results for Ekin + Eint , calculated with 
the exact ground state solution of the Gross-Pitaevskii 
equation, provides a test of the validity of the TF ap- 



proximation. For N ~ 10^ the agreement with the exact 
result is only semiquantitative, becoming better and bet- 
ter as N increases |^. After long expansion time the 
release energy can be related to the value of the square 
radii of the system, through the equation 
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The release energy has been measured both at Jila |T^] 
and MIT |||. It is worth noticing that in the MIT trap, 
where the initial configuration is a strongly anisotropic 
ellipsoid with the major axis along z, the release energy 
is almost entirely converted into kinetic energy of the 
radial motion, the velocity along z being much smaller. 
This can be tested by solving equation ( |l3| ) with the pa- 
rameters appropriate for the trap of Ref. Q . One finds 
that, after an expansion of 40 ms, the 2-body interaction 
energy is a factor 10^* smaller than the initial value and 
the ratio between the axial and radial kinetic energies is 
approximately 4 x 10"'^. 

The same formalism can be used to investigate the os- 
cillations of the trapped gas. One can easily check that, 
in the limit of small deformations, the solutions of ( |lO| ) 
yield the dispersion relation discussed in Ref. [|l^ . For an 
axially deformed trap the normal modes are classified in 
terms of the third component m of the angular momen- 
tum. We will discuss here the m = and m = 2 modes 
which are accounted for by the parametrization (|^^. 
The m — 2 mode, in the linear limit, has the frequency 
Lu = ^/2ujj_, while the low- lying m = mode, result- 
ing from the coupling between monopole and quadrupole 
oscillations, has the frequency ||l^ uj^ — rjuj'^, with 
7/ = (4-I-3A2-V9A4 - 16A2 + 16)/2. The frequency of the 
collective modes is expected to change when the ampli- 
tude of the oscillations becomes large, due to nonlinear 
effects. In order to calculate such deviations, we solve 
(|lo| ) using, as initial conditions, the ground state values 
bi{Q) — 1, but with 6(0) ^ 0. We choose the values of the 
velocities 6^(0) in order to to excite, in the linear limit, 
the two separate m — 2 and m — modes. For the m = 2 



mode this implies 6j. 



e and 6^ = 0, where the 



parameter e fixes the amplitude of the oscillations. For 
the m = mode, one has b^ = by = e and b^ = e{ri — 4). 
In this case, the system oscillates also along z, the axial 
and radial oscillations having relative amplitude ry — 4. It 
is worth noticing that the occurrence of a simultaneous 
oscillation in both the radial and axial widths is a typical 
effect of the interaction between the atoms. In fact, in the 
absence of 2-body forces (noninteracting harmonic oscil- 
lator), the motion in the two directions would be exactly 
decoupled. The experiments of Ref. |^ reveal not only a 
good agreement with the predicted frequencies of the two 
modes but also a clear evidence for the coupling between 
the axial and radial oscillations (see Fig. 2 of Ref. [^). 
These results are crucial signatures of the important role 
played by the interaction. 
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By increasing the initial values hi, it is possible to ex- 
plore the nonlinear regime. A major problem for the 
comparison with the experimental data is that the oscil- 
lations are measured by imaging the atomic cloud after 
switching off the trap and leaving the atoms to expand 
for a few ms. During the expansion the relative ampli- 
tude of the axial and radial motions can be significantly 
modified. This is especially true for the MIT trap, due to 
the strong asymmetry of the starting cigar-shaped con- 
figuration which makes the expansion in the radial and 
axial directions quite different. This nontrivial evolution 
of the oscillating cloud is not surprising if one thinks that, 
after an expansion time of a few ms, the size of the sys- 
tem can increase by more than a factor ten. In order to 
make a significant comparison with the experiments it is 
then important to simulate both the oscillations in the 
trap and the subsequent expansion. 

In Fig. ^ we show our predictions for the frequencies 
of the TO = and to = 2 modes, in the case of the Jila 
trap (A = VS), as a function of the relative amplitude, 
defined as {l/2)[{x^)UL- {x^)lll]/ {x^)lil. The latter is 
calculated after expanding the trap for 7 ms, as in the ex- 
periments of Refs. H,^. We find that the frequency of the 
TO = mode does not exhibit any significant dependence 
on the relative amplitude, while the to = 2 frequency in- 
creases. This agrees with the experimental findings [§,D, 
though the measured frequency shift of the m — 2 mode 
is about a factor two larger than our prediction. 

As already pointed out, the relative amplitude of the 
oscillations in the radial and axial directions behave dif- 
ferently during the expansion. For example, after excit- 
ing the TO = 2 mode in the trap, with a relative ampli- 
tude of 5% and 4% in the radial and axial directions re- 
spectively, the corresponding oscillations of the expanded 
cloud, after 7 ms, have a relative amplitude of 9% and 
4%. The aspect ratio Rr after 7 ms is found to be 1.65. 
These predictions rather well agree with the experiments 
(see Fig. 2 of Ref. |^), where the aspect ratio is found to 
be 1.75 and the radial and axial relative amplitudes are 
10% and 4%, respectively. 

We have repeated the same calculations for the cigar- 
shaped trap of the MIT experiments (A = 0.077 §1). In 
this case we excite the low-lying to = mode in the trap 
and then we calculate the frequency and the amplitude 
of the oscillations after expanding the gas for 40 ms. We 
do not observe any frequency shift of the low-lying m = 
mode as a function of the amplitude, in agreement with 
the experiments. One should note that the nonlinear- 
ity of the oscillations are amplified during the expansion 
much more than in the Jila trap, since the system is very 
anisotropic. As a consequence, the large axial ampli- 
tude observed after an expansion of 40 ms corresponds 
to a rather small amplitude at t — and this explain in 
part the absence of shift in the frequency. As an exam- 
ple, let us consider the aspect ratio plotted in Fig. 2 of 



Ref. [H, which oscillates between about 0.28 and 0.38. 
One can easily reproduce the same oscillations by solv- 
ing (|l^) starting from a configuration which oscillates as 
a low- lying m = mode in the trap (see also Ref. p6t). 
To get the measured aspect ratio one has to start with an 
oscillation having relative amplitude less than 1% in the 
radial direction and about 3.3% in the axial one. During 
the expansion the relative amplitude of the radial motion 
remains practically unchanged while the one of the axial 
motion increases up to 15 %. These large fluctuations 
of the axial width of the system are compatible with the 
conservation of energy, since almost all the energy is car- 
ried by the motion of the atoms in the radial expansion, 
which is much faster than the axial one. 

It is also worth noticing that, when A <C 1, the fre- 
quency of the high- lying to = mode is exactly 2ll!± even 
in the nonlinear regime. In this limit, this mode corre- 
sponds to a two-dimensional motion and the absence of a 
shift in its frequency reflects the occurrence of a hidden 
symmetry of a 2D Bose gas in a harmonic trap |2l[ . 

In conclusion we have shown that the hydrodynamic 
equations of superfluids provide a useful description of 
several nonlinear effects, associated with the dynamic be- 
havior of a trapped Bose gas at zero temperature. A rea- 
sonable agreement with the first available experimental 
data is found, though for a more quantitative compar- 
ison, when the number of atoms is relatively small, the 
complete solution of the Gross-Pitaevskii equations is ex- 
pected to be relevant. Our analysis points out the crucial 
role played by the interatomic forces in the dynamics of 
the expansion as well as in the behavior of the collective 
excitations. A natural extension of this work should in- 
clude thermal effects and, in particular, the interaction 
between the consensed and thermal components of these 
systems. 
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FIG. 1. Aspect ratio for the Jila trap (a;^ = 27rl32 Hz and 
A = \/8) as a function of the expansion time. Upper and lower 
solid lines correspond to 7?„ and Rr, respectively. Dashed 
lines are the same quantity for the noninteracting particles. 
Points are experimental values of Rr taken from pc(]. 



FIG. 2. Aspect ratio for the MIT trap {uji_ = 27r320 Hz 
and A = 0.056) as a function of the expansion time. Lower 
and upper solid lines correspond to f?„ and Rr, respectively. 
Dashed lines are the same quantity for the noninteracting 
particles. Points are experimental values of Rr extracted from 
Ref. [& (see text). 



FIG. 3. Asymptotic aspect ratio vs. A = uoz/l^_l- Solid 
lines: from (|l3|); dashed line: noninteracting gas. 



FIG. 4. Frequencies of the collective modes as a function 
of the relative radial amplitude, for the m = (triangles) and 
m = 2 (circles) modes in the Jila trap. Solid lines: from (^o|); 
points: from Refs. 
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